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Abstract 

The m a 6 correction to energy is expressed in terms of an effective Hamiltonian for an arbitrary 
state of helium. Numerical calculations are performed for n = 2 levels, and the previous result for the 2 3 P 
centroid is corrected. While the resulting theoretical predictions for the ionization energy are in moderate 
agreement with experimental values for 2 3 S±, 2 3 P, and 2 l So states, they are in significant disagreement for 
the singlet state 2 l P\. 
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High precision calculations of helium energy levels including relativistic and QED effects is 
a complicated task yj]. It has been recognized early on that the two-electron Dirac-Coulomb 
Hamiltonian is only an approximate Hamiltonian, as it includes negative energy spectra and does 
not account for magnetic and higher order interactions between electrons. The proper approach 
has to be based on Quantum Electrodynamic theory. For heavy few-electron ions the interactions 
between electrons can be treated perturbatively, on the same footing as the electron self-energy and 
vacuum polarization. Highly accurate results have been obtained for heavy helium- and lithium- 
like ions [2J, |3|], and a convenient formulation of this 1/Z expansion has been introduced a few 
years ago by Shabaev in |4|]. For systems with a larger number of electrons the zeroth order 
Hamiltonian will include an effective local potential to approximately account for interactions 
between electrons. This approach is being developed by Sapirstein et al |5|], and more recently 
by Shabaev and collaborators |(]]. One of the most interesting results obtained so far was the 
calculation of QED corrections to parity violation in the cesium atom [Q]. 

For light atomic systems relativistic and QED effects are only a small correction to the nonrel- 
ativistic Hamiltonian, and for this reason they can be treated perturbatively. More precisely, this 
perturbative approach relies on expansion of the binding energy in powers of the fine structure 
constant a 

E(a) = E® + E W + £ (5) + £ (6) + £ (7) + 0(a 8 ), (1) 

where E^ = ma n £^ is a contribution of order a n . However, this expansion is nonanalytic, 
inasmuch as some of the £^ coefficients contain In a, see for example Eq. ©. Each can be 
expressed in terms of the expectation value of some effective Hamiltonian H^ n > with the nonrel- 
ativistic wave function [|70. This approach allows for a consistent inclusion of all relativistic and 
QED effects order by order in a. We present in this work high precision calculations of n = 2 
energy levels in helium including the contribution E^ 6 ' . This contribution has already been derived 



1 111 . Here we obtain valid 



separately for triplet states in (8|, |90, and for singlet states in ifiol 
for all helium states, and present numerical results for 2 3 Si, 2 3 P, 2 1 Sq and 2 l P\ energy levels. 
The leading term in the expansion of the energy in powers of a, 8^ = £, is the nonrelativistic 



energy, the eigenvalue of the nonrelativistic Hamiltonian, which in atomic units is 
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The relativistic correction £ ^ is the expectation value of the Breit-Pauli Hamiltonian IU2I1 
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£( 5 ) is the leading QED correction. Apart from the anomalous magnetic moment correction to 
the spin-orbit and spin-spin interactions, which we neglect here, as we consider singlet or spin- 
orbit averaged (centroid) levels, it includes the following terms QJLXI 
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The next order contribution £^ is much more complicated. It can be represented in general as 
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but separate matrix elements of the first and the second term in the above are divergent. The spin 
dependent terms which contribute to fine structure are finite, and have been derived by Douglas 
and Kroll in [8|]. These contributions are not included here, because we consider spin-orbit av- 
eraged levels. The singularities of matrix elements in Eq. © can be eliminated by algebraic 
transformations [11] in a similar way for both singlet and triplet states. Therefore we extend the 
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result obtained in llllh to arbitrary states of helium, and the contribution £^ can be represented as 
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where P = pi + p 2 , p = (pi — p 2 ) /2, f ' = r\ — r 2 , and 
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The operators if^, ife, H c and ifo are parts of the Hamiltonian from Eq. ©, which was 
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transformed [11] to eliminate singularities from second order matrix elements 
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where V 2 V 2 in H' A in Eq. (flOl) is understood as a differentiation of on the right hand side as 
a function (omitting <5 3 (r)). and E'po are one- and two-loop electron self-energy and vacuum 
polarization corrections respectively [j 1 0L 1 1 If] - 
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The higher order contribution £ ^ is known only to some approximation. Following Ref. [13] 
the hydrogenic values for one-, two-, and three-loop contributions |l4l] at order ma 7 are extrapo 
lated to helium, according to 

(5 3 (r 1 ) + 5 3 (r 2 )) Hc 



£^ = [£( 7 )(LS,He + ) +^ (7) (nX,He H 
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for X = S, and for states with higher angular momenta £ ( - 7 \nX, He + ) is neglected. 

We pass now to the calculation of matrix elements. The wave function is expressed in terms of 
explicitly correlated exponential functions 0j 
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with random a*,/?*, 7* Ill5ll . This basis set is a very effective representation of the two-electron 
wave function, so much so that the nonrelativistic energies with 1500 basis functions are accurate 
to about 18 digits. Moreover, matrix elements of operators for relativistic and higher order correc- 
tions can all be obtained analytically in terms of rational, logarithmic and dilogarithmic functions, 
for example 

d 3 r l / d 3 r 2 = — -. (19) 



16 7r 2 J J ' r 1 r 2 r (a + (3){f3 + + a) 

Numerical results for matrix elements of m a 6 operators with singlet and triplet P states are pre- 
sented in Table I: due to the singularity of these operators we had to use octuple precision arith- 
metic. The m a 6 correction to the energy also involves second order matrix elements E sec . If we 
write 

H' A = Q A (20) 
H B = Q B -s (21) 
Hc = Qc- { ^pt (22) 

H D = Q%s l s> (23) 

then one obtains for singlet states 

E(2 1 S ) sec = (2 1 S\Q A Q A \2 l S) + (2 1 3\Q 3 C Q^S) (24) 

E(2 1 P 1 ) SCC = (2 1 P i \Q A g A |2 1 P i ) + (2 1 P i |^ c Q^P 1 ) (25) 

and the contributions from Hb and Hd vanish. The result for the 2 3 5i state is 

E(2 3 S 1 ) SCC = (2 3 S\Q AJ ^Q A \2 3 S) + ^(2 3 S\Q BJ ^Q B \2 3 S) 

+ I (2 3 S\Q 3 C Q J C \2 3 S) + I (2 3 S\Q% — L_ Q%\2 3 S) (26) 

The result for the 2 3 P centroid, defined by 

E{2 3 P) = l [E(2 3 P ) + 3 E{2 3 P 1 ) + 5 E{2 3 P 2 )] (27) 
9 

is 

E(2 3 P) sec = (2 3 P*\Q Al ^Q A \2 3 P l ) + ^(2 3 P l \Q B ^^Q B \2 3 P<) 

+ I (2 3 P l \Q J c J^e Q ^ PI) + I ^ P ^d Qg|2 3 P*> (28) 



TABLE I: Expectation values of operators entering for the 2 l P\ and 2 3 P centroid. 
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-0.039 975 


0.071 600 




0.973 055 


1.198 492 




3.102 248 


3.883 404 


pi x p 2 /rpi x p 2 


0.216 869 


0.399 306 


p^ p^ ^ 1 r f'^ 7^ j f"~* | 3 j ^5 ^ pi pJ 


-0.126416 


-0.187304 
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TABLE II: Contributions to ionization energ y £ ( 6 > for n = 2 states of the helium atom. En is a sum of 
operators in Eq. <[Sj> , in comparison to Ref. jllll it includes the contribution Eg. E^q is the logarithmic 
contribution, last term in Eq. ([8]). The sum of spin dependent, second order corrections for 2 3 P centroid is 



taken from Ref. llSfl. 



m a 


2 X S 


2 X P 


2 3 S 


2 3 P 


E Q 


12.287491 


12.236 966 


13.052 109 


11.963 305 


E' A 


-16.280186(10) 


-16.084034(5) 


-17.189 809(10) 


-15.848 510(2) 


Eb 


0.0 


0.0 


-0.018 722 . 




Ec 


-0.033 790 


0.201363 


-0.001 108 


> - 0.168 704(2) 


Ed 


0.0 


0.0 


-0.003 848 ' 




-P D irac(He + ) 


4.000000 


4.000000 


4.000000 


4.000 000 


Subtotal 


-0.026485(10) 


0.354 296(5) 


-0.161377(10) 


-0.053 908(3) 


Em 


2.999960 


0.106839 


3.625 397 


-1.106 416 


Er2 


0.016860 


0.000112 


0.032 331 


-0.009 867 


Elg 


0.133682 


0.011364 


0.0 


0.0 


Total 


3.124017(10) 


0.472 611(5) 


3.496 351(10) 


-1.170191(3) 



Numerical results for second order matrix elements are presented in Table II. One notices a strong 
cancellation between m a 6 contributions and the Dirac energy for the He + ion, the subtotal line in 
Table II. Because of this cancellation, the dominant contribution is the one loop radiative correc- 
tion, with the exception of the 2 1 P 1 state, where the wave function at the nucleus happens to be 
very close to 16, the He + value, see Table I. 

The summary of all important contributions to ionization energies is presented in Table III. We 
include the first and second order mass polarization correction to the nonrelativistic energy, as well 
as first order nuclear recoil corrections a 4 m 2 /M and a 5 m 2 /M. We expect higher order terms in 
the mass ratio to be much below the 0.01 MHz level, the precision of calculated contributions, see 
Table III. Results for nonrelativistic as well as for leading relativistic corrections are in agreement 



with those obtained previously by Drake 1113 



0,fl. 



Corrections of order ma 5 were calculated 
using the Drake and Goldman Ill7ll values for Bethe logarithms. The m a 6 correction is calculated 
in this work. All but m a 7 contributions are calculated exactly. This last one, m a 7 is estimated 
on the basis of the hydrogenic value according to Eq. (fl~6l) . It is the only source of uncertainty 
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TABLE III: Contributions to ionization energy of n = 2 helium states in MHz. Physical constants from [23], 
Roo = 10 973 731.568 525(73) m -1 , a = 1/137.03599911(46), ^ e = 386.1592678(26) fm, m a /m e = 
7294.2995363(32), c = 299792458. E is is a finite nuclear size correction with the charge radius r a = 1.673 
fm. 











uC? 3 P \ 


2 

fia - 


-960 331428.61 


-814 736 669.94 


-1152 795 881.77 


-876 058183.13 




8 570.43 


41522.20 


6711.19 


-58 230.36 


u 3 IM 2 a 2 


— 16 72 


—20 80 


— 7 1 1 


—25 33 


Efs 


1.99 


0.06 


2.59 


-0.79 


A 

m a 


-11971.45 


-14 024.05 


-57629.31 


11436.88 


m 2 /Ma 4 


-3.34 


-2.81 


4.28 


11.05 


m a 5 


2 755.76 


38.77 


3 999.43 


-1234.73 


m 2 /Ma 5 


-0.63 


0.47 


-0.80 


-0.62 


m a 6 


58.29 


8.82 


65.24 


-21.83 


m a 


-3.85(1.90) 


-0.16(16) 


-5.31(1.00) 


1.93(40) 


Ethe ~ 


-960332 038.13(1.90) 


-814 709147.44(16) 


-1152 842 741.56(1.00) 


-876106 246.93(40) 




-960332 040.86(15) 


-814 709153.0(3.0) 


-1152 842 742.97(0.06) 


-876106 247.35(6) 



of theoretical predictions, as the achieved numerical precision for each correction is below 0.01 
MHz. 



The value for the 2 1 S'o state has already been presented in our former work f| 1 1|1 - here we 
display in more detail all the contributions. The value for the 2 3 Si state is in agreement with our 
previous calculation in B9D, where we obtained S^ 1 = 3.496 93(50). This provides justification 
of the correctness of the obtained result, since the two derivations of the m a 6 operators were 
performed in a different way. However, the result for the 2 3 P state is in disagreement with our 
result from [ISO. For this reason we checked Ref. [|18|l . and found a mistake. The derived set of 
operators representing 8^ was correct, but the expectation value of H' EN , in the notation of [18], 
was in error. The correct result is (H' EN ) = 11.903 751. With the second order matrix element 
-15.838 656(9) and subtracting He+ m a 6 energy -Z 6 /16, it is equal to 0.049 702(9), while the 
former result was 0.140 689(9), see Table II of [18]. Together with other corrections from that 
Table the total ma 6 contribution becomes —1.170 188(9), in agreement with the result obtained 
here. 

We find a moderate agreement with experimental ionization energies for the 2 1 S , 2 3 Si and 2 3 P 



states but a significant disagreement for the 2 1 Pi state. Following 111311 . the result for the 2 3 Si state 
was obtained by combining the 2 3 Si -3 3 L>i measurement by Dorrer et al UM 786 823 850.002(56) 
MHz with the theoretical 3 3 D 1 ionization energy 366 018 892.97(2) MHz calculated by Drake 
ifl^ flfjl. The ionization energy of the 2 3 P state was obtained from the measurement of the 2 3 S' 1 — 
2 3 P transition by Cancio et al [20] of 276 736 495.6246(24) MHz and the previously obtained 
2 3 Si ionization energy. The ionization energy of the 2 1 S' state was obtained from measurements 
of 2 1 S — n l D transitions by Lichten et al 112111 with n — 7 — 20 and Drake's calculations for 
l D states QQ]. Finally, the result for 2 1 P ionization energy is determined by combining the 
2 1 P-3 1 P 2 transition 448 791 404.0(30) MHz by Sansonetti and Martin D22D (including correction 



of 0.6 MHz ITJ), with calculated 113, 



]M 3 1 J D 2 energy 365 917 749.02(2) MHz. The significant 



disagreement with theoretical predictions for 2 1 P state calls for an independent calculation of 
the m a 6 term, and on the other hand for the direct frequency measurement of 2 1 P — 3 1 D 2 or 
2 1 P - 2 ^transitions. 

Further improvement of theoretical predictions can be achieved by the calculation of ma 7 
contributions. The principal problem here will be the numerical evaluation of the relativistic cor- 
rections to Bethe-logarithms and the derivation of remaining operators. Such a calculation has 
recently been performed for helium fine structure Il24ll . therefore in view of newly proposed ex- 
periments Il25ll . calculations for other states of helium although not simple, can be achieved. 
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